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Abstract 

In this work we show that the physics reach of a long-basehne (LBL) neutrino oscil- 
lation experiment based on a superbeam and a megaton water Cherenkov detector can 
be significantly increased if the LBL data are combined with data from atmospheric 
neutrinos (ATM) provided by the same detector. ATM data are sensitive to the octant 
of ^23 and to the type of the neutrino mass hierarchy, mainly through three-flavor effects 
in e-like events. This allows to resolve the so-called 622,- and sign(Am|]^)-parameter 
degeneracies in LBL data. As a consequence it becomes possible to distinguish the nor- 
mal from the inverted neutrino mass ordering at 2a CL from a combined LBL+ATM 
analysis if sin^ 2^13 > 0.02. The potential to identify the true values of sin^ 2^13 and 
the CP-phase 5cp is significantly increased through the lifting of the degeneracies. 
These claims are supported by a detailed simulation of the T2K (phase II) LBL exper- 
iment combined with a full three-flavor analysis of ATM data in the HyperKamiokande 
detector. 
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1 Introduction 

Thanks to the out-standing developments in recent neutrino physics a rather clear picture of 
neutrino oscillation parameters is emerging. We know that there are two large angles in the 
lepton mixing matrix, 612 and ^^23! and one small angle ^13, and the mass-squared differences 
Am^i and |Am|i| are well determined by the global data from solar [1] atmospheric [2], reac- 
tor [3] , and accelerator [4] neutrino experiments. A recent global analysis of world neutrino 
oscillation data can be found e.g. in Ref. [5]. Despite these enormous achievements there are 
still important open questions related to neutrino physics, which have to be addressed by 
future neutrino oscillation experiments: 

1. What is the value of the mixing angle Oi^l 

2. What is the value of the complex phase (5cp in the mixing matrix? 

3. Is the neutrino mass ordering normal (mi < m^. < m^) or inverted {m^ < mi < 777.2), 
ie., what is the sign of Amg^? 

In contrast to the "discovery phase", which was dominated by natural neutrino sources 
such as the sun or the atmosphere, the subsequent generation of oscillation experiments 
will be mainly based on man-made neutrinos, where the neutrino source is well under con- 
trol. Among the future projects are conventional beam experiments [6], new reactor exper- 
iments [7], supcrbeam experiments [8,9], and eventually experiments based on a neutrino 
factory or a beta-beam [10]. 

A characteristic feature in the analysis of future long-baseline (LBL) experiments is the 
presence of so-called parameter degeneracies, see e.g.., Refs. [11-17]. Due to the inherent 
three-flavor structure of the oscillation probabilities, for a given experiment in general several 
disconnected regions in the multi-dimensional space of oscillation parameters will be present. 
Traditionally these degeneracies are referred to in the following way: 

• The intrinsic or (5cp; 6*13) -degeneracy [18, 19]: For a measurement based on the — >• 

oscillation probability for neutrinos and anti-neutrinos two disconnected solutions 
appear in the (5cp, ^13) plane. 

• The hierarchy or sign(Am3j^)-degeneracy [20]: The two solutions corresponding to the 
two signs of Amg]^ appear in general at different values of 5cp and ^13. 

• The octant or ^^23-degeneracy [21]: Since LBL experiments are sensitive mainly to 
sin^26'23 it is difficult to distinguish the two octants ^^23 < t^/'^ and 6'23 > 7r/4. Again, 
the solutions corresponding to ^23 and 7r/2 — ^23 appear in general at different values 
of 5cp and ^13. 

This leads to an eight-fold ambiguity in the determination of 9i^, and 5cp [11]i and hence 
degeneracies provide a serious limitation in the ability to answer the questions 1 and 2 
above. Moreover, the fact that one speaks of a "sign(Am3^)-degeneracy" illustrates that 
answering question 3 is difficult: determining the neutrino mass ordering is equivalent to 
resolving the sign(Am|i)-degeneracy. Several methods to resolve these degeneracies have 
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been proposed, among them the combination of experiments at various basehnes and/or 
(L/i?) -values [11,20,22-25], the use of spectral information [12,26], the combination of 
and z/g oscillation channels [27], and the combination of LBL and reactor 

experiments [28-32]. 

In the present work we discuss a new possibility to resolve the LBL parameter degenera- 
cies, based on the data from atmospheric neutrinos. It is known that atmospheric neutrinos 
are in principle sensitive to ^13 and the neutrino mass hierarchy due to earth matter ef- 
fects [33] in the e-like events [34-40]. These effects are most pronounced in the multi-GeV 
energy range and for large zenith angles, corresponding to neutrino trajectories crossing the 
earth mantle or the mantle and the core. In addition effects from the solar parameters Qyi 
and Am2i on e-like events in the sub-GeV energy range provide sensitivity to the octant of 
023 [41-43], and in principle even on 6cp [44]. For a recent discussion of sub- leading effects 
in atmospheric neutrino oscillations see Ref. [45]. It turns out, however, that atmospheric 
(ATM) data on its own can never compete with LBL in many respects, such as the determi- 
nation of 6'i3,5cp, |^"^3il; or sin^2^23- One reason is that they are limited by systematical 
uncertainties, e.g. from the neutrino fluxes. However, as we will show in the following, due 
to the effects mentioned above ATM data can break the sign(Am3j^)- and ^23- degeneracies in 
LBL data, and hence the combined analysis of LBL and ATM leads to significant synergies. 

An important part of the future neutrino program are gigantic water Cherenkov detectors 
at the megaton scale. Apart from serving as detector for LBL experiments such a facility 
will provide unprecedented opportunities for proton decay, solar and atmospheric neutrino 
experiments, as well as for the detection of supernova and other astrophysical neutrinos. The 
projects under discussion are UNO [46] in the US, a megaton detector in the Frejus labora- 
tory [47] in Europe, and the HyperKamiokande project [9,40] in Japan. If a LBL experiment 
with such a detector will be built atmospheric neutrino data come for free. Therefore, our 
method provides a very efficient possibility to resolve the parameter degeneracies from LBL 
data, in contrast to the previously discussed methods, which in general are based on the 
combination of two or more expensive experiments. 

In the following we will illustrate how the LBL+ATM combination works by considering 
the phase II of the T2K experiment (T2K-II) [9], assuming a 4 MW superbeam produced 
at the J- PARC accelerator, and the 1 Mt HyperKamiokande (HK) detector serving as the 
far detector for the LBL experiment as well as providing the high statistics atmospheric 
neutrino data. Let us note that similar results are expected also for the other megaton 
detector proposals. Furthermore, if huge magnetized iron detectors are available one expects 
similar synergies between LBL and ATM data. For these type of detectors the ability to 
distinguish between neutrinos and anti-neutrinos will increase the potential of ATM data 
even further [48]. In the present work we stick to the HK water Cherenkov detector, the 
potential of magnetized detectors will be considered elsewhere. 

The outline of the paper is as follows. In Section |21 we give some technical details of our 
simulation of the T2K-II experiment f Section 12.11) and of the atmospheric neutrino analysis 
(Section l2.2|l . and we describe the statistical methods used in the following (Section l2.3|l . In 
SectionElwe illustrate the LBL/ATM complementarity by discussing the effects of parameter 
degeneracies for the LBL data f Section 13.11) . whereas in Sections 13.21 and 13.31 we consider the 
effects in ATM data which allow to resolve the octant and the sign(Am3;^) degeneracies, 
respectively. In Section HI we investigate in detail the ability to exclude degenerate solutions 
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The T2K-II LBL experiment 



The HK ATM experiment (9 Mt yrs) 



V (2 Mt yrs) 



V (6 Mt yrs) 



V 



Ve signal 
Uf. background 
signal 

background 



21300 
2140 
73 200 
340 



16 000 
3 260 
75 600 
320 



e-like sub-GeV 
e-like multi-GeV 
/i-like sub-GeV 
/i-like multi-GeV 
upward going // 



239 000 
52 700 
232 000 
108 000 
127000 



58 000 
18 100 
66 200 
49 100 
65 400 



Table 1: Number of events in the LBL and ATM experiments considered in our work for the osciUation 
parameters sin^ 26*13 = 0.05, sin^ 6*23 = 0.5, sin^ 6112 = 0.3, 5cp = 0, Am^^ = 8.1 x 10"^ eV^, and Am^;^ = 
2.2 X 10-3 eV2. 



with LBL+ATM data by performing a systematic scan of the parameter space of ^13, 6'23 
and ^cp- In Sectional we discuss the potential to identify the neutrino mass ordering, for 
the special case 6^23 = 7r/4 f Section 15. 1|) as well as for the general case f Section l5.2|) . In 
Section ini we show how the sensitivity to ^13 is improved by resolving the ^23-degeneracy by 
ATM data, and we summarize our results in Section [7| 

2 Description of the experiments and statistical analysis 
2.1 The T2K-II long-baseline experiment 

In the following the label "LBL" refers to the phase II of the T2K experiment (T2K-I I) [9]. We 
are assuming a high luminosity superbeam with mean neutrino energy of 0.76 GeV, produced 
with a target power of 4 MW at the J-PARC accelerator. The neutrinos are detected at a 
1 Mt water Cherenkov detector, HyperKamiokande (HK), at a baseline of 295 km and an 
off- axis angle of 2°. We consider 2 years running time with a neutrino beam, and 6 years 
with anti-neutrinos, such that comparable event numbers are obtained for neutrinos and 
anti-neutrinos. 

Our simulation of the T2K-II experiment is performed by using the GLoBES software 
package [49]. We take into account realistic neutrino fluxes, detection cross sections, energy 
resolution, and efficiencies [9]. This experimental information is folded with the three- flavor 
oscillation probability, fully taking into account the earth matter effect. We consider all 
available information, from —>■ v^. appearance, as well as z/^ disappearance channels. The 
signal events are given by z/g and charged current (CC) interactions, respectively. We 
divide the signal into the total event rate, where the full CC cross section is used, and into 
the energy spectrum with free normalization, where only events from quasi-elastic scattering 
are used, since non-quasi-elastic events do not allow to reconstruct the neutrino energy. 
Various backgrounds such contamination of the beam, neutral current events, 

and misidentification of muon neutrinos as electron neutrinos are taken into account, using 
information given in Ref . [9] . We list the numbers of signal and background events expected 
in T2K-II for typical oscillation parameters in Tabled More details of our T2K-II analysis 
can be found in Ref. [13]. 
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2.2 The HK atmospheric neutrino experiment 

To simulate atmospheric neutrino data in the HK detector we follow closely Ref. [43]. A 
crucial element of the analysis is to take into account the full three-flavor oscillation prob- 
ability, including earth matter effects, as well as oscillations induced by the "solar" mass 
splitting Amgi (for other three- flavor analyses see Refs. [50-53]). In our analysis we consider 
charged current data, divided into sub-GeV and multi-GeV e-like and //-like contained event 
samples (each grouped into 10 bins in zenith angle), as well as stopping (5 angular bins) and 
through-going (10 angular bins) up-going muon events. 

Details of our statistical analysis can be found in the Appendix of Ref. [54]. Together 
with the statistical errors, we consider theoretical as well as systematical uncertainties, where 
theoretical uncertainties are uncertainties in the original atmospheric neutrino fluxes and in 
the cross-sections. We are using the atmospheric neutrino fluxes from Ref. [55], and flux un- 
certainties include total normalization errors (20%) allowing for an independent fluctuation 
of neutrino and anti- neutrino fluxes (5%) as well as z/^ and fluxes (5%), spectral uncer- 
tainty of the fluxes ("tilt" factor), and an uncertainty on the zenith angle dependence which 
induces an error in the up/down asymmetry of events (5%). We also include independent 
normalization errors for the different contributions to the interaction cross section: quasi- 
elastic scattering (15%), single pion production (15%), and deep inelastic scattering (15% for 
contained events and 10% for upward-going muons). Moreover, we include as systematical 
errors experimental uncertainties associated with the simulation of the hadronic interactions, 
the particle identification procedure, the ring-counting procedure, the fiducial volume deter- 
mination, the energy calibration, the relative normalization between partially-contained and 
fully-contained events, the track reconstruction of up-going muons, the detection efficiency 
of up-going muons, and the stopping/through-going separation.^ 

The current atmospheric neutrino data sample from SuperKamiokande (SK-I) consists 
of 1489 days of data (contained events) with a detector mass of 22.5 kt, which gives roughly 
90 kt yrs of data. In this work we are considering a data taking period for the LBL experiment 
of 8 years. We assume that the 1 Mt HK detector will be finished one year before the T2K-II 
beam, and hence 9 Mt yrs of atmospheric neutrino data will be available. For the atmospheric 
data sample used in our analysis, which we will denote in the following by the label "ATM" , 
we scale the present SK-I sample (1489 days contained events, 1657 days of stopping, and 
1678 days of through-going muons) by a factor 100. Event numbers for various ATM data 
samples for typical oscillation parameters are given in Table 

2.3 Details of the statistical analysis 

In order to investigate the potential of the experiments described in the previous sub-sections 
we adopt the standard method for analyzing future experiments. First, artificial "data" is 
simulated by calculating event numbers for LBL and ATM for some "true values" 0*''™ for 
the oscillation parameters 6 = (^13, ^12, 6'23, 5cp, Am2]^, Amg^). Then a x^-analysis of these 
"data" is performed to extract allowed regions for the oscillation parameters. It is important 
to note that in general the results will depend on the values adopted for 0^^^°. For all 

^The impact of these theoretical and systematical uncertainties on the performance of future atmospheric 
neutrino experiments has been investigated in some detail in Ref. [43] . See also Ref. [45] . 



5 



calculations we will use for (Am2i)*'^"'^ and (Am|^)*''"^ the best fit values obtained in 
Ref. [5], 

sin2^i2 = 0.3, Am^^ = 8.1 X 10-^eV^ jAm^^] = 2.2 x IQ-^eV^ (1) 

We do not expect any significant changes of our results if these parameters are varied within 
the present allowed ranges [5]. However, we adopt various values for Of^^jO^^'^, 5^™^, and 
sign(Am3^), and show our results as a function of these parameters. Furthermore, since 
neither LBL nor ATM data allow an accurate determination of 6*12 and Am^i we assume 
that these two parameters are known with an uncertainty of 10% from solar and reactor 
neutrino experiments. 

If one is interested in the allowed range for a certain parameter ^ G for a given choice 
of O^^'^^, the function x^i^^^'^^] ^) has to be minimized with respect to all other parameters 
d except ^ to take into account the correlations and degeneracies between parameters. This 
minimization is performed by using the GLoBES software package [49], which has been 
generalized in order to include the atmospheric neutrino code. Let us stress that for both 
data samples, LBL as well as ATM, a full three-flavor analysis including matter effects is 
performed. The only approximation is to neglect the (very weak) dependence on 9i2. Since 
LBL depends in leading order only on the product sin 26*12 Am^i fixing 6*12 does not introduce 
any error, as long as the dependence on Am2i is properly taken into account, and for ATM 
varying 9 12 in the allowed range is expected to have a very small impact. Apart from this 
simplification the full parameter dependence of both data sets has been taken into account. 
However, we assume that LBL and ATM data are statistically independent. This might 
not be completely correct since both data sets are based on the same detector, and hence, 
uncertainties related to the detection process can introduce correlations between LBL and 
ATM data. Although such effects will have to be included eventually in the analysis of real 
data, we expect that such correlations introduce only very minor corrections to our present 
results. 

3 Complementarity between LBL and ATM data 

3.1 Parameter degeneracies and the T2K-II experiment 

In this section we discuss the problem of parameter degeneracies for the T2K-II experiment. 
First we note that for this particular experiment the (5cp, ^i3)-degeneracy does not occur. 
It is known that for experiments operating at the oscillation maximum of Amg^ the sec- 
ond solution in the (5cp,6'i3) plane can be disfavored efficiently [11,13]. Moreover, spectral 
information is important, since (5cp-dependent and 5cp-independent terms in the oscilla- 
tion probability show a different energy dependence (compare Eq. Q later in this section). 
Thus it is difficult to leave the spectrum unchanged when 5cv is varied, and the (5cp,6'i3)- 
degenerate solution, which implies a different value of (5cp, is disfavored. An illustration of 
the relevance of spectral information for the (5cp; 6'i3)-degeneracy in the T2K-II experiment 
can be found in Figure 5 of Ref. [13]. 

Therefore, in the case of T2K-II we are confronted only with the sign(Am|]^)- and ^23- 
degeneracies. Apart from the incapacity to determine the neutrino mass ordering and the 
octant of 6*23 this leads to a four-fold ambiguity in the determination of 6*13 and (5cp. The 
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True sinV =0.4 True sin\ = 0.6 




Figure 1: Allowed regions in the (sin^ 26*13, (5cp) plane at 2cr, 99%, and 3a CL (2 dof) of the true and 
all degenerate solutions for sin^ 2ef^'' = 0.03, S^^j^'' = -O.SStt, and sin^ fl^^j™ = 0.4 (left) and sin^ 9^^^'= = 0.6 
(right). The solid curves correspond to LBL data only, and the shaded regions correspond to LBL+ATM 
data. The true best fit point is marked with a star, the best fit points of the degenerate solutions are marked 
with dots, and the corresponding Ax^-values of LBL+ATM data are given in the figure. The true mass 
ordering is the normal hierarchy. 

impact of tlie degeneracies can be appreciated in Figure ^ where the solid curves show the 
allowed regions from LBL data in the (sin^ 2^13, 5cp) plane for an example-point with the 
true values sin^ 2^13 = 0.03, ^cp = — O.SSvr, and non-maximal values of O^^'^. Apart from 
the true solution, three degenerate regions are present, corresponding to the wrong octant 
of 6*23, the wrong sign of Amg^, and the wrong octant as well as the wrong hierarchy. Let us 
introduce the following abbreviations to denote these four solutions: 

true solution 

true octant of 623 and wrong hierarchy , . 

wrong octant of ^23 and true hierarchy 
wrong octant of 623 and wrong hierarchy 

A qualitative understanding of the degenerate solutions can be obtained from the ap- 
proximate formula for the z/^ —>■ Ve appearance probability in vacuum [24] 

P,^_,^ ~ sin^ 2^13 sin^ ^23 sin' A 

+ a sin 26*13 sin 26*12 sin 26^23 A sin A cos( A ± (5cp) (3) 
+ cos' ^23 sin' 2^12 A' 

with A = Aml^L / (AEu) and a = Am2i/Am|i. The sign in the second term is '+' for 
neutrinos and ' — ' for anti-neutrinos. Since in the T2K-II experiment matter effects are small 
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this expression for the probabihty suffices to obtain a quahtative understanding of most of 
the effects in LBL data presented throughout this work. 

Following Ref. [11] the location of the 6'23-degenerate solution (O™'', if*'') can be estimated 
from Eq. (jS)) in the following way. Since the T2K-II experiment operates close to the oscil- 
lation maximum it is a good approximation to use A ^ 7r/2. Furthermore, for values of 
sin^25i3 > 0.01 the term can be neglected. Under these approximation solving the two 
equations 

Pu.^.Ms, Sep, 623) = P.,^.,(0i3,(5^p,7r/2-^23) 

Pp,-.P.(^13,5CP,^23) = Pp,^p.(^;3,^CP>^/2-^23) (4) 

leads to 

sin^ 20^3 ^ sin^ 2^13 tan^ 623 (5) 

sin(5cp ^ sin^cp tan 6^23 (6) 

for the parameter values of the wrong-6'23 solution. For the example shown in Figure C] 
Eq. dSD gives sm^2e[^ = 0.02 for sin'^ 6%'''' = 0.4, and sin2 2^'i3 = 0.045 for sin^ 6%'''' = 
0.6. Furthermore, Eq. © gives S'^p = -O.Slvr for sin^ O^^f^ = 0.4, and ^cp = O.SSvr for 
sin^ 6'23'^'^ = 0.6. These numbers are in good agreement with the actual fit shown in Figure ^ 
One observes that the ^23- degeneracy has a strong impact on the measurement of sin^26'i3, 
however, it hardly affects the determination of 6cp [11,13]; the (O™'', ii*'^)-solution occurs 
practically at the same value of 6cp as the true one. 

This is in contrast to the sign(Am3^)-degeneracy, which in general leads to a severe 
ambiguity for 6cp, whereas the sin^ 2^13 measurement is essentially unaffected. Following 
Ref. [20] we observe from Eq. (jS)) that only the term in the second line is affected by changing 
the sign of Amg^ . This term transforms as 

a A sin A cos(A ± 6cp) ^ -a A sin A cos(-A ± 6cp) (7) 

under Amg^^ —>■ — Amg;^. Therefore, the full probability stays invariant under this transfor- 
mation if 6cp is adjusted such that cos(— A ± ^cp) = — cos(A ± Sq^), which can be achieved 
for 

(5^p = vr - Sqp . (8) 

According to this reasoning we would obtain in the example plotted in Figure ^ the value 
6'qp = 1.857r= — O.IStt. Although in this case the accuracy is not excellent, Eq. (jHJ still 
provides a rough method to estimate the location of the (0*^^, ii™'')-solution in ^cp (for more 
refined methods see e.g. Refs. [14,56]). Finally, the values of sin^ 2^13 and 6cp corresponding 
to the combined 6*23 and sign(Am|j^)-degeneracy {0^^,H'"^) can be estimated by applying 
simultaneously Eqs. (0) and 

This four-fold degeneracy can be lifted to large extent if LBL data is combined with 
data from atmospheric neutrinos. We observe from Figure ^ that the degenerate solutions 
corresponding to the wrong octant of ^23 are highly disfavored by the inclusion of ATM 
data, at the level of Ax^ ^ 20. Furthermore, also the solution with the wrong mass ordering 
gets disfavored in the combined analysis, although in this case the ability to resolve the 
degeneracy is more subtle. In the following subsections we discuss the relevant features of 
ATM data to resolve the degeneracies in LBL data. 



8 



3.2 The sensitivity of ATM data to the octant of 6*23 



The T2K-II experiment will provide a very precise determination of sin^2^23 thanks to the 
large statistics data from the i/^ disappearance channel (compare Table ^). The relative 
accuracy at la is expected to be better than 1%. Despite this impressive performance on 
sin^ 26^23 there are some subtleties related to the measurement of sin^ ^23 (see, e.g.^ Refs. [57, 
58]). Especially, if ^23 deviates from 7r/4 it will be impossible to distinguish between the two 
solutions at ^23 and 7r/2 — ^^23-^ Although the determination of sin^26'23 from atmospheric 
data is significantly less precise than from LBL they provide the very interesting ability to 
distinguish between the two octants of ^23- 

The sensitivity of atmospheric data to the deviation of 6'23 from 7r/4 follows mainly from 
effects of the solar mass splitting Am2x for e-like events in the sub-GeV region, see e.g. 
Refs. [41-44]. The excess of e-like events can be written as 

= ^-1 ^ (r COS^ ^23-1) (Pin . (9) 
e 

Here Ag (Ag ) is the number of e-like events with (without) oscillations, {P2\) is the averaged 
two neutrino oscillation probability given by the solar parameters Amgi and Qyi including 
the (weighted) contributions from neutrinos and anti-neutrinos, and r = F^/F^ is the ratio 
of the initial muon and electron neutrino fluxes. Since for sub-GeV energies r ^ 2 the effect 
is suppressed for ^23 ~ 7r/4, however it provides a sensitive measure for deviations from 
maximal mixing. Given the present LMA parameters, for | sin^ ^23 — 0.5| ^0.1 one expects 
gsub yg^i^gg Qf a, few percent (see e.g. Fig. 8 of Ref. [44]). In contrast to LBL data, which 
is essentially sensitive only to sin^2^23; this effect depends on cos^^^23; and therefore the 
discrimination between 6'23 > 7r/4 and ^23 < 7r/4 becomes possible. 

Building upon the results of Ref. [43] we show in Figure El the difference in between 
the true solution (with = 0) and the x^-minimum in the wrong octant of ^23 ((O™'', H^'^)- 
solution). It is clear form this flgure that LBL data alone have no sensitivity at all to the 
octant of ^^23, whereas ATM can reject the fake solution efliciently. Taking the combination 
of LBL+ATM data improves the sensitivity slightly, due to the more precise determination 
of other oscillation parameters by the LBL data. Using the combined LBL+ATM data the 
wrong octant can be rejected at 3a if | sin^ 6*23 — 0.51 > 0.1. Let us note that the results shown 
in Figure El do hardly depend on our choice for the true neutrino mass ordering, thanks to 
^13"*^ = 0. Non-zero values for 6\'^^^ will be considered in Section |31 

3.3 The sensitivity of ATM data to the mass hierarchy 

The determination of the ordering of the neutrino mass states is one of the most challenging 
tasks of future neutrino experiments. In long-baseline experiments the effects of 6*13, 6cp 
and the sign of Amgj^ are highly entangled, and the determination of sign(Am3^) is probably 
only possible through the matter effect induced in an experiment with a very long baseline 
plus additional information from other LBL or reactor experiments [13, 20, 22, 23, 30, 60- 

^For theoretical expectations for the deviations of 623 from the maximal value see e.g. Refs. [57, 59] and 
references therein. 
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Figure 2: Ax^ of the solution in the wrong octant of ^23 as a function of the true value of sin^ 623 for 
LBL data only, for ATM data only, and for the LBL+ATM combination. Furthermore, we take ^^3"° = 0. 



62].'^ Therefore, the combined analysis of LBL and ATM data as offered from an oscillation 
experiment with a Mt water Cherenkov detector provides a very interesting possibility to 
answer the question of the neutrino mass hierarchy. 

The sensitivity of atmospheric neutrino data to the neutrino mass hierarchy comes mainly 
from the modification of e-like multi-GeV events by earth matter effects for not too small 
values of sin^ 2^13 [34-40]. Similar as in Eq. 0, one finds for the 6*13 induced excess of e-like 
multi-GeV events 

6r'" = ^ - 1 ^ (r sin^ - 1) (Pir) . (10) 

e 

Now (P^i) is an effective two-flavor probability governed by Am|]^ and ^13, appropriately 
averaged and including the (weighted) contributions from neutrinos and anti-neutrinos. The 
effect is most pronounced for zenith angles corresponding to neutrino trajectories crossing the 
earth mantle, or earth mantle and core, where sin^ 26'13-effects can be resonantly enhanced 
due to matter effects [33-35]. In the relevant zenith angle bins g™*^^*' can reach values of 
the order of 10% (see e.g. Fig. 5 of Ref. [39]). Quahtatively 6™"i*' shows the following 
behavior [39]: 

• vanishes for sin^ 2^13 = and increases monotonically with sin^ 2^13. 

• For the normal hierarchy the resonant matter enhancement occurs for neutrinos, whereas 
for the inverted hierarchy it occurs for anti-neutrinos. Since the event numbers in water 
Cherenkov detectors are dominated by neutrinos because of larger cross sections (see 
Table P), e™"'*i is larger by a factor of 1.5 — 2 for the normal hierarchy than for the 
inverted one. 

■^Another possibility to identify the neutrino mass hierarchy could come from the observation of neutrinos 
emitted by a galactic supernova, see e.g. Ref. [63] for a recent analysis. 
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• Since in the multi-GeV range we have r ~ 2.6—4.5, the factor (r sin^ ^'23 — 1) in Eq. (|TT1|) 

suppresses the excess of e-hke events for small values sin^^23 ^ 0.4, whereas larger 
values of sin^ ^23 increase e™"!*'. In particular, the effect is enhanced for ^23 > 7r/4. 

Therefore, if the true hierarchy is normal and ^23 < 7r/4 there is only a small excess of e-like 
multi-GeV events, which can be accommodated to some extent with the inverted hierarchy. 
For the example point with sin^ 9^^'^ = 0.4 chosen in the left panel of Figure ^ the wrong 
hierarchy can be disfavored only with a Ax^ = 5.0. In contrast, for 6'23 > vr/4 the excess of 
e-like events is enhanced by the flux-factor (r sin^ ^23 — 1)- Therefore, a true normal hierarchy 
(resonant enhancement for neutrinos) plus the fiux-factor enhancement leads to large values 
of e™"'*', which cannot be fitted with the inverted hierarchy (resonant enhancement for the 
smaller anti- neutrino sample). This explains the strong rejection of the wrong hierarchy 
solution in the right panel of Figure ^ with a Ax^ = 18.6. We will discuss the mechanisms 
relevant for the rejection of the wrong hierarchy in more detail in Sections lU and El 

Finally we note that in addition to Am^^-effects of Eq. ^ and 6'13-effects of Eq. (fTO|l 
also an interference term between the two contributions is present [44] (see also Ref. [53]). 
It is proportional to (r sin ^13 sin 26*23) and depends on the CP-phase ^cp- Because of the 
different dependence on the flux ratio r the interference term may become important in cases 
where the effects governed by Eqs. and (fTIH) are suppressed. 

4 Resolving the degeneracies 

The examples shown in Figure suggest that the LBL-I-ATM combination offers an efficient 
method to reject the degenerate solutions. To investigate this more systematically we have 
performed a scan over the true values for ^13 and 6*23 and various values for Sqp'^. For a 
given point in the space of sin^ 2^*™^, sin^ ^23"'' and ^^p*^ ^^st the ability to rule out each 
of the three degenerate solutions (O^^if*''), {O^^H""'), and (O™^ //"■•)• This is done by 
minimizing Xlbl+atm with respect to all fit parameters, constraining the octant of 6*23 and 
the sign of Amli corresponding to the fake solution which we want to test. The results of 
this analysis are shown in Figure El which represents one of the main results of this work. 

The solid curves in Figure El show that the (O™'', if*'')-solution corresponding to the ^23- 
degeneracy can be excluded at high confidence level if ^23 is far enough from 7r/4. (Note 
that if ^23 is close to 7r/4 this degeneracy disappears anyway.) For small values of sin^ 29^^"^ 
this follows mainly from the atmospheric sub-GeV e-like events, as discussed in Section 13.21 
Moreover, if the true hierarchy is normal (left panels of Figure E]) we find an improvement 
of the octant sensitivity for sin^ 2^13 > 0.04. This effect comes from the multi-GeV e-like 
events, where resonant enhancement occurs for sufficiently large sin^ 2^13, and therefore an 
additional dependence on sin^6'23 appears according to Eq. (fTIUl . This can be seen from 
Figure m where we show the ability to reject the fake solutions using only sub-GeV (left) 
and multi-GeV (right) ATM data. The improvement for large sin^ 2^13 is also visible for 
the inverted mass ordering (right panels of Figure EI); however, in that case the effect is 
much smaller, since the resonance occurs for anti-neutrinos in the inverted hierarchy, which 
contribute less to the total ATM data, and hence the significance of the effect is smaller than 
in the normal hierarchy, where the resonance occurs for neutrinos. 

The wrong hierarchy solution {O^^, H""^) can be excluded to the right of the dashed curves 
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Figure 3: Rejection of the fake solutions from LBL+ATM data as a function of the true values of sin^ 29i3 
and sin^ 623 for various Sq™ if the true hierarchy is normal (left) or inverted (right). Solid curves correspond 
to the solution with the wrong octant and the right hierarchy (O^'', H^'^), dashed curves to the right octant 
and the wrong hierarchy (O*'', 77'^''), and shaded regions to the wrong octant and the wrong hierarchy 
(O™'^ , H™'^). We show the contours of Ax^ — 1,4,9 between the fake and the true solution, corresponding 
to a rejection of the fake solution at the la, 2a, and 3a CL (from light to dark shading) for 1 dof. 
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Figure 4: Rejection of the fake solutions from LBL data combined with sub-GeV (left) and multi-GeV 
(right) e-like and /i-like ATM data as a function of the true values of sin'^ 26*13 and sin^ 6*23 for d^p"^ — 
(upper) and (5^™° ~ — 7r/2 (lower) and true normal hierarchy. Solid curves correspond to the solution with the 
wrong octant and the right hierarchy (O'"'', i?*"^), dashed curves to the right octant and the wrong hierarchy 
(O'"^, 7?^''), and shaded regions to the wrong octant and the wrong hierarchy {O^'^ , H^'^). We show the 
contours of Ax^ — 1,4,9 between the fake and the true solution, corresponding to a rejection of the fake 
solution at the la, 2a, and 3a CL (from light to dark shading) for 1 dof. 



in Figure El We observe that the abihty to reject this solution increases with sin^023- As 
discussed in Section 1331 this follows from the fact that the excess of e-like multi-GeV events 
is enhanced for large sin^^23- In particular, for sin^ ^23 > 0.5 the sensitivity is completely 
dominated by multi-GeV data (compare Figure Ej), and we observe a very small dependence 
on the true value of 6cp and on the true hierarchy. For sin^ ^23 < 0.5 the situation becomes 
more complicated, and the sensitivity depends on the true value of 6cp- In that region also 
sub-GeV data can contribute significantly, as visible in FigurelH For low values of sin^ 6^23 the 
effect of multi-GeV data becomes suppressed by the flux factor (r sin^ ^'23 — 1) in Eq. (fTIUl . and 
the contribution of sub-GeV data from Eq. can become of comparable size. Moreover, 
as discussed in Ref. [44] for certain parameter values an interference term between the two 
contributions may become relevant, which depends on 6cp- Therefore, the final sensitivity 
emerges from a rather involved interplay of various effects and ATM data samples. 

Finally, moving to the (O™"^, if'"'')-solution, we find that in most cases this "double" 
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degenerate solution can be excluded if either the (O™"^, H^"^)- or the (O*'', if'^'')-solution is 
ruled out. However, for certain parameter configurations an interesting degeneracy between 
the octant and the mass ordering appears, such that the (O*'', if™'') -solution cannot be 
excluded, although the (O™'', H^^)- and the (O*'', if™'^)-solutions are not allowed. We find 
from Figure El that it is difficult to distinguish a normal mass hierarchy and sin^ 6'23 — 0.45, 
5cp ~ from an inverted hierarchy and sin^6'23 — 0.55, (5cp — ^r. We observe from the 
upper panels of Figure Elthat for these particular parameters sub-GeV as well as multi-GeV 
data cannot exclude the (O™'', if™'')-solution. However, the appearance of this remaining 
degeneracy strongly depends on the true value of ^cp- For example, from the lower panels 
of Figure m it becomes clear that it does not occur for 5cp = — vr/2. The dependence on 5cy> 
indicates that again the interference term between Am2i- and 6'13-effects is important. 

5 Determining the neutrino mass hierarchy 

In this section we discuss in detail the possibility to identify the type of the neutrino mass 
ordering from a combined LBL+ATM analysis. To this aim we simulate data for the T2K-II 
LBL experiment as well as for atmospheric neutrino data assuming a "true" neutrino mass 
ordering. Then we fit these data with the "wrong" hierarchy and search for the minimum 
X^-value in the full 6-dimensional space of oscillation parameters, taking into account both 
^23-degenerate solutions (O*'', H""^) and (O™"^, H""^). If this minimum is larger than a certain 
value the wrong hierarchy can be excluded at the corresponding CL. 

5.1 Sensitivity to the mass hierarchy for maximal 6*23 mixing 

Before considering the general case we will first adopt the choice 9^2z'^ = 7r/4, such that 
the 6'23-degeneracy is absent, and we are left only with the two- fold ambiguity related to 
sign(Am3;^). We show the sensitivity to the hierarchy in Figure El as a function of the true 
values of sin^ 26'i3 and 5cp. First we observe from this figure that the sensitivity of LBL 
data alone strongly depends on the true value of 6cp- For certain values of 5cp the ability 
to distinguish normal and inverted hierarchy is even lost for sin^ 26'i3 = 0.1. The main 
reason for the difficulties of LBL alone to determine the mass hierarchy, is that for T2K-II 
the matter effect is very small. ATM data on their own allow to identify the normal mass 
hierarchy at 2a CL for sin^ 2^^™*^ > 0.04 (see left panel of Figure E)), whereas there is rather 
poor sensitivity to the inverted hierarchy (right panel of Figure jSl). However, the sensitivity 
is significantly increased by combining the two data sets: For LBL+ATM data the wrong 
hierarchy can be excluded at 2cr CL for sin^ 2^*3"^^ > 0.02, with a rather small dependence 
on the true value of 6cp or the true type of the hierarchy. 

These results can be understood qualitatively from the discussion given in Section 13.31 
To achieve the same value of e™"'*i as implied by the (true) normal hierarchy and a given 
sin^ 2^^3™ with the inverted hierarchy, one has to adopt values of sin^ 2^13 larger than 
sin^ 29^^^ to increase the effect for the inverted hierarchy. Conversely, if the true hierar- 
chy is inverted and one wants to fit e™"^*' with the normal ordering, sin^ 2^13 has to be 
smaller than sin^ 2^*3^'°. This expectation is confirmed in Figure IHl where the allowed regions 
in the (sin^ 2^13, 6cp) plane from ATM data are shown as dashed curves for the true values 
sin^ 2^13 = 0.03 and 6cp = — 0.857r. The left panel shows that the normal hierarchy can be 



14 




Figure 5: Sensitivity to the mass hierarchy as a function of the true values of sin^ 26'i3 and Sqp for 
023^° = 7r/4, if the true hierarchy is normal (left) or inverted (right). We show the contours of Ax^ = 1,4,9 
between the wrong and the true hierarchy, corresponding to a rejection of the wrong hierarchy at the la, 
2a, and 3a CL (from light to dark shading) for 1 dof. The shaded regions correspond to LBL+ATM data 
combined, solid curves correspond to LBL-only, and dashed curves to ATM-only. 



fitted with tlie inverted one for rather large values of sin^ 26'i3, whereas the opposite situation 
is visible for the inverted hierarchy in the right panel. 

This behavior explains also the ATM-only results shown in Figure El The reasonable 
sensitivity of ATM data to the normal hierarchy visible in the left panel is given by the fact, 
that a fit with the inverted hierarchy would be only possible for very large values of sin^ 2^13. 
This is disfavored by data, first because the values of sin^ 2^13 required to fit the data start 
to get in confiict with the CHOOZ bound, and hence are disfavored, and second because 
such large values of sin^ 2^13 start to be excluded from ATM data on their own. On the 
other hand the reason for the poor sensitivity to the inverted hierarchy shown in the right 
panel of Figure |S1 is that in general small values of e™"'*^ are expected, which can be easily 
fitted with the normal hierarchy but smaller sin^ 2^13 (compare right panel of Figure IHl), well 
below the sensitivity of ATM data. 

Let us now discuss the effect of the sign(Am3]^)-degeneracy for LBL data. In agreement 
with the discussion of Section ITTl we observe from Figure IHl that the sign{Am1i) degeneracy 
affects mainly the determination of 6cp, whereas the value of sin^ 2^13 of the wrong-sign 
solution practically coincides with the true one. According to Eq. (jH)) we expect for the fake 
solution the value Sqp = — O.IStt, which is in good agreement with the values obtained in 
the actual fit: 6qp = — 0.25(— 0.05)7r for true hierarchy normal (inverted). The deviations 
from the value of 6qp given by Eq. (jH)) are due to the matter effect, which is not included in 
the probability Eq. (jH)) and to some extent also to spectral information, which in general is 
difficult to include in the analytical discussion. Note that because of the choice 6*2™^ = 7r/4 
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Figure 6: Allowed regions in the (sin^ 26'i3, i^cp) plane assuming for the true mass spectrum a normal (left) 
and an inverted (right) hierarchy, and 6*2™° = 7r/4. The true values sin^ 26*53 



0.03 and (5^^^^= 



-O.SStt 



are marked by a star. The shaded areas and the solid curves correspond to the 1, 2, and 3 cr allowed regions 
(2 dof) for the true and the sign(Am32)-degenerate solution for LBL+ATM and LBL data, respectively. The 
dashed curves correspond to the Ax^ contours of ATM data for the wrong sign of Ato§i . The Ax^ of the 
best fit point with the wrong hierarchy for LBL+ATM data is given in the figure and its location is marked 
by a dot. 



the (O*'', H""^)- and (O™'', if"'')-solutions coincide, and we are left only with the two-fold 
ambiguity implied by the sign(Am|i)-degeneracy, as visible in Figure IHl 

From Figure IHl also the reason for the significant improvement of the sensitivity to the 
hierarchy for the LBL+ATM combination becomes clear. Although for LBL data alone the 
degenerate solution is present at a wrong value of Sqp even at la CL, the value of sin^ 2^13 
has to be very similar to the true one. This makes it impossible to accommodate ATM data 
with the wrong hierarchy, since either values sin^ 29i3 significantly larger or smaller than the 
true one are necessary to obtain the correct value for e^'ii*' with the wrong hierarchy. The 
location in the (sin^ 26^13, 5cp) plane required by the degenerate solution of LBL gives for 
ATM data a corresponding to 2 — 3cr CL. Therefore the combined data lead to the rather 
good sensitivity implied by the Ax^-values given in Figure IHl 

To summarize, the example shown in Figure IHl demonstrates how the combination of LBL 
and ATM data leads to a sensitivity to the mass hierarchy at the 3a level, although each 
data set on its own cannot distinguish normal from inverted mass ordering for this particular 
choice of parameter values. The sensitivity of ATM data alone suffers from the fact that 
the wrong hierarchy can be accommodated by adjusting 613, but also the parameters 6*23, 
Amg^, and ^cp are important. To benefit from the earth matter effects in e-like events from 
atmospheric neutrinos for the hierarchy determination the rather precise measurement of the 
oscillation parameters from LBL data is mandatory. In particular, it is necessary to know 
023 and Am|^ at the 1% level, and a reasonable constraint on sin^ 2^13 is required. Finally, 



16 




Figure 7: Sensitivity to the mass hierarchy as a function of the true values of sin^ 2^13 and sin^ for 
5q^° — and a true normal (left panel) and inverted (right panel) hierarchy. We show the contours of 
Ax^ — 1, 4, 9 between the wrong and the true hierarchy, corresponding to a rejection of the wrong hierarchy 
at the Icr, 2(7, and 3(t CL (from light to dark shading) for 1 dof. The shaded regions correspond to LBL+ATM 
data combined, solid curves correspond to LBL-only, and dashed curves to ATM-only. 



the determination of tlie correct mass hierarcliy through the LBL+ATM combination allows 
to resolve the ambiguity in 5cp from the sign(Am3j^)-degeneracy, and hence the potential to 
measure (5cp is significantly increased. For the example shown in Figure El the correct value 
of (5cp can be identified at 2a CL. 



5.2 Sensitivity to the mass hierarchy for non-maximal ^^23 mixing 

Now we generalize the discussion of the previous subsection to non- maximal values of ^23"*', 
such that in general all four solutions of Eq. (0) are present. To identify the correct mass 
hierarchy the (O*'', H""^)- as well as the (O™'', if™'')-solution has to be excluded. In Figure [7| 
we show the sensitivity to the mass hierarchy as a function of the true values of sin^ 2612, 
and sin^ ^^23 for LBL, ATM, and LBL+ATM data. First we note that the results for LBL 
data-only strongly depend on the true value of 5cP) as observed already in Figure ATM 
data on their own provide only a reasonable sensitivity to the mass hierarchy if the true 
hierarchy is normal and sin^ 6*23'^'^ > 0.5. This is the region of large excess of e-like multi-GeV 
events (see discussion in Section EIHl) , which cannot be achieved for any configuration in the 
inverted hierarchy. 

The sensitivity to the hierarchy for the LBL+ATM combination can be inferred from 
Figure El by considering the intersection of the shaded region ((O™"^, iJ™'^)-solution), with 
the region to the right of the dashed curves ((O*'', i7'"'')-solution). The shaded regions in 
Figure [71 give an explicit example for ^^p'^ = 0. For the true inverted hierarchy (right 
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Figure 8: Panel (a) shows the Ax^ of the wrong hierarchy, panel (b) and (c) show the values of sin^ 2^13 
and sin^ ^23 at the best fit point of the wrong hierarchy as a function of the true value of sin^ ^23- The true 
hierarchy is normal and sin^ 26^;^° = 0.04, S^^^^" = 0. 



panel) the sensitivity to determine the mass hierarchy is given by the abihty to exclude the 
{O^^, i?'^'") -solution. As discussed in Section El the worsening of the sensitivity for low values 
of sin^ can be understood from the flux factor (r sin^ ^23 ~ 1) in the multi-GeV e-like event 
excess. If the true hierarchy is normal (left panel) the (O*'', iJ"'')- as well as the {O^^, H""^)- 
solutions have to be taken into account. The "spike" visible in Figure [3 at sin^ 6'23"^ ~ 0.45 
comes from the fact that the double degenerate solution (O™'', if"'') cannot be excluded. 

We illustrate this behavior in Figure |Hl where we show the location of the best fit point 
for the wrong hierarchy as a function of the true value of sin^ 623- We observe from panel (c) 
that for sin^ ^23^'^ > 0.5 the minimum for LBL+ATM follows the true value of sin^^23; ^-e., 
it corresponds to the (O*'', if'"'')-solution. As discussed in Section El this implies in turn 
that sin^ 2^13 has to be close to the true value, as can be seen from panel (b). However, 
in that case the large values of e™"'*i from the normal hierarchy cannot be obtained in the 
inverted hierarchy, which leads to the large values, i.e. good sensitivity, in that region. For 
0.4 < sm^923^ < 0.5 one observes from Figure |H1 that the best fit moves to the [O^'^.H^^)- 
solution, characterized by the wrong ^23? see panel (c), and a value of sin^ 26'i3 given by 
Eq. (0), see panel (b). As mentioned in Section HI the occurrence of this degeneracy strongly 
depends on the true value of (5cp, and follows from a dehcate interplay of effects in sub- 
and multi-GeV data. In the region of sin^ ^2™° < 0-4 the best fit returns again to the 
(O*'', if'^'')-solution, characterized by the true 622, and a value of sin^ 26*13 close to the true 
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Figure 9: Sensitivity to the mass hierarchy as a function of the true values of sin^ 26*13 and Sqp for 
sin^ ^2™° — 0-3 (left) and sin^ ^23"° — 0-7 (right). The true hierarchy is normal. We show the contours of 
Ax^ — 1, 4, 9 between the wrong and the true hierarchy, corresponding to a rejection of the wrong hierarchy 
at the la, 2a, and 3a CL (from light to dark shading). The shaded regions correspond to LBL+ATM data 
combined, solid curves correspond to LBL-only, and dashed curves to ATM-only. 

one. From Figure |H1 we find tliat tlie fit of LBL-only, as well as ATM-only chooses best 
fit values of 6*23 and sin^ 26*13 close to the true ones. This suggests that in that region the 
sensitivity comes from effects related to Amg^^, Arrigi, and 6cp- As mentioned in Section |3] 
here the main sensitivity comes from the ATM sub-GeV data sample (compare Figure 0} and 
the interference term between Am^^- and Amg^-effects [44] can become important, which 
introduces the dependence on ^cp- 

The strong dependence of the sensitivity on 6cp for small sin^ 6*23 is visible in the left panel 
of Figure El where we show the sensitivity to the hierarchy as a function of sin^ 291^^^ and 
Sqp' for sin^ ^23^'^ = 0.3. The rich structure visible for the LBL-I-ATM combination indicates 
a complicated interplay and/or cancellations of various effects. On the other hand, from 
the right panel we observe that for large values of sin^ ^23"'^ the sensitivity to the hierarchy 
becomes practically independent of the true value of 6cp- In this range multi-GeV Amg^- 
effects are enhanced by the flux factor and dominate over Am2i-effects and the interference 
term. Also ATM data alone provide quite a good sensitivity, thanks to the big effect in 
multi-GeV e-like events. 

6 The ^23-degeneracy and the sensitivity to sin^ 2^13 

As discussed in Sectional the presence of the octant degeneracy leads to severe complications 
in the determination of sin^ 26*13. In this section we consider the specific case ^13"*^ = 0, 
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Figure 10: Sensitivity to sin^ 2^13 as a function of the true value of sin^ ^23 for LBL data only (dashed), 
and the combination LBL+ATM (solid). In panels (a) and (b) we restrict the fit of 6'23 to the octant 
corresponding to 623^'' and 7t/2 — 62'^^°, respectively, whereas panel (c) shows the overall sensitivity taking 
into account both octants. 



and we investigate how the upper bound on sin^ 2^13 ("sensitivity to sin^2^i3") can be 
improved by the combined LBL+ATM analysis. The fact that the 6'23-degeneracy affects 
the sin^ 26'13-sensitivity is well-known, see e.g., Refs. [11,13,14,16]. The choice ^^3"*^ = 
essentially eliminates the differences between normal and inverted hierarchies, i.e. the 
degeneracy is perfect, and the fact that it cannot be resolved introduces only very small 
ambiguities for other parameters. Hence the (O™'', if*'')- and (O""^, if'"'')-solutions coincide 
within good accuracy. Moreover, the analysis is independent of the true value of 6cp, since 
the phase becomes unphysical for 6*13 = 0. We define the sensitivity to sin^ 26*13 as the largest 
value of sin^ 26*13 which fits the data generated for a true value sin^ 26**3"*^ = at a given CL 
(see Appendix C of Ref. [31] for a detailed discussion). 

Let us first assume that the true octant of 6*23 was known. In that case one expects from 
Eq. (jni) that the sensitivity for sin^ 2^13 becomes better for increasing sin^ 623, since for large 
sin^ 6*23 the first term in Eq. (jSJ, which is proportional to sin^2^i3, is enhanced. Moreover, 
the term gets suppressed for small cos^ 623 , which reduces the effect of multi-parameter 
correlations induced by Am^i. This expectation is confirmed in panel (a) of Figure^l where 
the sensitivity to sin^ 2^13 is shown as a function of the true value of sin^ 6*23, assuming that 
the octant is known. However, if the octant of 623 is not known, the degenerate solution 
prevents the smaller sin^ 2^13-limits implied by sin^ 6*2™^ > 0.5. In that case the data can be 
fitted by relatively large values of sin^ 26*13, since exchanging sin^ 623 and cos^ 6*23 reduces the 
effect of sin^ 26*13 and increases the term in Eq. (jS)). This effect is shown by the curves 
corresponding to LBL data in panel (c) of Figure ^1 We observe that the sensitivity to 
sin^ 2^13 gets worse for non-maximal values of 6*23"*^ in both octants because of the presence 
of the degenerate solution. 

Since the octant- degeneracy can be efficiently resolved by atmospheric data one expects a 
significant improvement of the sin^ 2^13-sensitivity for the combined LBL+ATM analysis. In 
panel (b) of Figure^! we show the sin^ 2^13-limit constraining the fit to the wrong octant of 
6*23. As expected from Figure |21 there are no solutions for LBL+ATM data if 6*23 is sufficiently 
far from maximal, since the wrong solution is disfavored by ATM data. Consequently, we 
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find that the final result shown in panel (c) for LBL+ATM is very close to the situation in 
panel (a), where only the true solution for 623 has been used. Hence, the sin^ 2^13-sensitivity 

is significantly increased for 623 > 7r/4 through the exclusion of the wrong octant solution by 
ATM data. There is no relevant improvement for 623 < ir/A, since in that case the sensitivity 
to sin^ 29i3 suffers for the true 6*23 solution because of the sin^ 26'13-suppression implied by 
the small values of sin^ 623. 

7 Conclusions 

In this work we have performed a combined analysis of future long-baseline (LBL) and atmo- 
spheric (ATM) neutrino data. As a specific example we have considered the phase II of the 
T2K experiment, consisting of a 4 MW superbeam produced at the J-PARC facility. The 
1 Mt HypcrKamiokande water Cherenkov detector will serve as detector for the LBL exper- 
iment and simultaneously provide high statistics ATM neutrino data. We have shown that 
the combined LBL-I-ATM analysis offers a very appealing possibility to resolve parameter 
degeneracies in the LBL data. In particular, the ambiguities implied by the sign(Am|i)- and 
the ^23-degeneracies can be lifted to a large extent. This becomes possible through three- 
fiavor effects in ATM data related to ^^13 and Am'^i. A systematic scan of the parameter 
space has been performed to investigate the ability to determine the type of the neutrino 
mass ordering. Let us summarize our main findings: 

• For true values of sin^ ^^23 > 0.5 the correct mass hierarchy can be identified at 2o" CL 
if sin^ 26*13 ^ 0.015, rather independent of the true value of 5cp or the type of the true 
hierarchy. In this region the sensitivity is dominated by multi-GeV e-like events in 
ATM data. 

• For true values of sin^ ^23 < 0.5 the correct mass hierarchy can be identified at 2a CL if 
sin^ 2^13 > 0.03, where here the actual sensitivity depends on the true value of 5cp and 
the true hierarchy. The final sensitivity emerges from an interplay of various effects in 
the different ATM data samples. 

• The solution with the wrong octant of ^23 can be excluded at Sa CL if | sin^ ^23 — 0.5| > 
0.1, independent of the true values of ^13, Scp, and the hierarchy. This follows mainly 
from Am2i-effects in e-like sub-GeV ATM data. If sin^ 2^^i3 > 0.03 and the true 
hierarchy is normal the octant sensitivity at 3a CL improves to | sin^ 6*23 — 0.5 1 > 0.05 
due to ^13-effects in the multi-GeV ATM data. 

• The lifting of the degeneracies by ATM data significantly increases the performance 
of the LBL experiment for the measurement of sin^ 26*13 and ^cp, since fake solutions 
implied by the degeneracies can be ruled out. Generically, the determination of the 
correct octant of ^^23 removes an ambiguity in the measurement of sin^2^i3, whereas 
lifting the sign (Amg^^) -degeneracy allows the identification of the correct value of Sqp- 

• If no finite value for ^13 is found the upper limit on sin^ 2^13 is significantly improved 
by resolving the ^23-degeneracy by ATM data if ^23 > 7r/4. 

Let us stress that these results follow from the complementarity of the two data sets; neither 
LBL data alone nor ATM data alone can provide a comparable physics reach. The LBL data 
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allow a very precise determination of |Am|^| and sin^2^23, and although ^13 and 6cp suffer 
from ambiguities related to the degeneracies they are constrained to rather specific values. 
However, there is very poor sensitivity to the mass hierarchy and to the octant of 623 for 
LBL data alone. For ATM data alone only for sin^ 62;^'^ ^ 0.6 and a true normal hierarchy 
a reasonable sensitivity exists for the mass ordering, since in this case a big excess of e-like 
multi-GeV events is predicted, which cannot be achieved by any configuration within the 
inverted hierarchy. For all other regions in the parameter space the sensitivity of ATM- 
only is rather poor, since the data can be fitted with the wrong hierarchy by adjusting the 
oscillation parameters. To benefit from the hierarchy sensitivity offered by ^13-earth matter 
effects in ATM data the precise measurement of the oscillation parameters from LBL data 
is mandatory. In particular, it is necessary to know 623 and |Am|^| at the 1% level, and a 
reasonable constraint on sin^ 2^13 is required. Hence, only the combination of LBL and ATM 
data makes it possible to obtain a good sensitivity to the neutrino mass hierarchy. Let us 
add that to benefit from the LBL+ATM combination indeed a few Mt yrs of ATM data are 
necessary. In particular we have checked that for exposures below 1 Mt yrs the sensitivity to 
the hierarchy is essentially lost. 

Finally, let us remark that the methods to resolve parameter degeneracies discussed pre- 
viously in general involve the combination of two or more (expensive) experiments, e.g. at 
different baselines or using different oscillation channels. In contrast, once a LBL experiment 
with a megaton Cherenkov detector is built, ATM data come for free. Therefore, in addi- 
tion to the interesting physics, the synergies between LBL and ATM data offer a relatively 
economical way of resolving parameter degeneracies. 
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